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ABSTRACT 



We present a pilot study of the stellar populations of 8 AGN hosts at z ~ 1 and compare to 
(1) lower redshift samples and (2) a sample of nonactive galaxies of similar redshift. We utilize 
K' images in the GOODS South field obtained with the laser guide star adaptive optics (LGSAO) 
system at Keck Observatory. We combine this K' data with B, V, i, and z imaging from the 
Advanced Camera for Surveys (ACS) on HST to give multi-color photometry at a matched spatial 
resolution better than 100 mas in all bands. The hosts harbor active galactic nuclei (AGN) as 
inferred from their high X-ray luminosities {Lx > lO**^ ergs s~^) or mid-IR colors. We find a 
correlation between the presence of younger stellar populations and the strength of the AGN, as 
measured with [OIII] line luminosity or X-ray (2-10 keV) luminosity. This finding is consistent 
with similar studies at lower redshift. Of the three Type II galaxies, two are disk galaxies and 
one is of irregular type, while in the Type I sample there only one disk-like source and four 
sources with smooth, elliptical/spheroidal morphologies. In addition, the mid-IR SEDs of the 
strong Type II AGN indicate that they are excited to LIRG (Luminous InfraRcd Galaxy) status 
via galactic starbursting, while the strong Type I AGN are excited to LIRG status via hot dust 
surrounding the central AGN. This supports the notion that the obscured nature of Type II AGN 
at z ~ 1 is connected with global starbursting and that they may be extincted by kpc-scale dusty 
features that are byproducts of this starbursting. 

Subject headings: galaxies: active (galaxies:) quasars: general galaxies: Seyfert instrumentation: 
adaptive optics techniques: high angular resolution X-rays: galaxies 



INTRODUCTION 



It is observed that the evolution of active galactic nuclei ( AGN) and quasi-stellar objects (QSOs) is 
connected with the propert ies of the bulges in their host galaxies ( Magorrian et al.lll998l : Ferrarese fc Merritt 
2OOOI : ICebhardt et al.ll2000l ) at low redshifts. These correlations between black hole mass and various bulge 
properties (e.g., velocity dispersion, half-hght radius, total mass) hint at simultaneous evolution of the AGN 
and its host galaxy. While this connection between activity at sub-parsec and kiloparsec scales is observable 
in snapshots of an AGN's life history, a complete understanding of its evolution through time remains elusive. 
It is also unclear how the condit ions and morphology o f the galaxy at large may be related to the canonical 
dust torus on sub-parsec scales ( Antonucci et al. 1993h . 



Several scenarios explaining these observations have arisen in recent years from simulations of galaxy 
formation. In the major merger scenario, two colliding gas-rich galaxies trigger central, dusty star formation 
and black hole fueling, followed by an unobscured quasar (QSO) phase resulting as AGN feedback expels 
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the gas and dust (e.g., Sanders et al. IQSSt Hopkins et al. 2005 : Springel et al. 2005h . This quenching of 
star formation may be a mechanism by which AGN activity controls the evolu tion of galaxies from the 
quiescent blue cloud to the red sequ ence jHopkins et all bood: ICroton et aPbood ). serving as a component 



of the "mixed" evolution scenario of 



Faberetal 



(|2007f) . An alternate picture is that the triggering of AGN 
activity simply depends on the conditions at the time. Both the condition of existence of a supermassive 
black hole (SMBH), met in most if not all ellipt icals, and the con dition of pres ence of a large reser voir of 



cold gas, are required for AGN activity to occur ( Kauffmann et al. 12003. 2007; S ilverman et al. 20081 ) 



Large, carefully controlled studies of the morphological properties, emission properties, and stellar 
populations of AGN hosts at low and high redshift have shed light on these theories. Several studies 
have found that luminous quasars and radio galaxies inhabit massive ellipticals formed at or before z > 2.5 
( McLure et al iNolan et al .I2OO1I ). Contrastingly, other studies have found that intermediate-luminosity 



AGN (1042 < Ljf,o.5- 



lOfcey 



< 10 ergs s ) at low redshift tend to reside in a unique class of objects that 



possess elliptical morphologies, but additionally display the signatures of young stellar po pulations (YSPs 



which places them in the green valley region between the blue cloud and red sequence. iKauffmann et al 



([20031) find that narrow-line AGN (Type II, 0.02 < z < 0.3) of all luminosities reside almost exclusively in 
massive galaxies and have distributions of sizes, stellar surface mass densities and concentrations that are 
similar to those of ordinary early-type galaxies. In this Sloan Digital Sky Survey (SDSS) study, the hosts of 
high-luminosity AGN have much younger mean stellar ages than those of the low-luminosity AGN, which are 
similar in stellar makeup to old ellipticals. In addition, the young stars are not preferentially located near the 
nucleus of the galaxy, but are spread out over scales of at least several kiloparsecs. The interpretation put 
forth is that these "young bulges" simply possess the necessary conditions for AGN activity - a SMBH and a 
reservoir of cold gas. Many other studies have uncqyered evidence of young s t ellar populat ions in l ocal Type 



II AGN JSchmitt et al]jl999t Icid Fernandes et all boOll Ijoguet et aP I2OO1I : iGonzalez Del gado etliD I2OOI 



Storchi-Bergmann et al.l 120011 ) . Th e tendency of AGN to reside in hosts with e lliptical-like morphologies 



has been seen in previous studies (jStorchi-Bergmann et al.l I2OOII : iHo et al. _ 2003 ). Add i tiona l information 



on the spatial distribution of YSPs in local AGN hosts is presented in IKauffmann et al.l (j2007l ). Hosts with 
young bulges and strongly accreting black holes often display outer disks (r > 10 kpc) with blue near-UV - 
r colors, suggesting that outer star formation is connected with AGN growth, according to near-ultraviolet 
(UV) fluxes from the Galaxy Evolution Explorer (GALEX) . The authors hypothesize again that these AGNs 
possess reservoirs of cold gas, th e raw ma terial for star formatio n. Galaxy interactions (Gun^Esi^ Hernquist 



1989; 



Hernquist fc Mihogll995l ) or bars (jShlosman et al.lll992l ) must then be invoked to funnel this gas into 



the nucleus and thence onto the black hole. 

Although the quantity of work on these issues at low redshift is substantial, AGN evolution theories 
additionally make predictions about the properties of their progenitors at higher redshift (z 1 — 2). In 
particular, it was once thought that the major merger - AGN fueling scenario predicted that mergers and 
close compaiiions w ould be more common in AGN hosts, but this has not been observed ([Pierce et al.ll2007 : 
Grogin et al ] |2005L Direct observations of AGN feedback have been difficult, as the timescale for this process 
is short (< 1 Gvr. [Silverman et al.ll2008i ). cutting down on number stati stics. S e veral studies have focused on 
the morphological properties of high redshift AGN hosts. For example, iHasanl ([2007^ finds that the hosts of 
intermediate luminosity AGN (10^^ < < 10"^^-^ ergs s""'^) at z ~ 1 are m ore likely to be massive ellipticals 



or bulge-dominated disks, just as at low-redshift. [Silverman et al.l (|2008r ) analyzes the rest-frame colors of 
the hosts of intermediate luminosity AGN at z ~ 1 and arrives at a similar conclusion. 

Studies of the stellar populations of AGN hosts at high redshift can impart independent information 
about their individual histories than can be inferred from their structural or morphological properties alone. 
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For example, the presence of YSPs in high-redshift AGN hosts would signal that the processes controlling 
AGN fueling were not vastly different from today. A comparison of the ages of dominant stellar populations 
for a sample of AGN against a sample of nonactive systems could indicate a delay time between AGN activity 
and starbursting. 

In this paper we search z 1 AGN hosts for young stellar populations, using adaptive optics imaging 
together with data from the Hubble Space Telescope Advanced Camera for Surveys (HST ACS) to enable 
stellar population modeling via 5-band Spectral Energy Distribution (SED) fitting. We investigate AGN in 
the Great Observatories Origins Deep Survey South (GOODS-South) field, tak ing advantage of the large 



number of AGN revealed by the 1 Megasecond Chandra exposure in this field (IGiacconi et al. 2001). We 



use the Laser Guide Star Adaptive Optics (LGSAO) system on the Keck Telescope (jWizinowich et al.ll2006l ) 
with the Near Infrared Camera-2 (NIRC2) at K'-band and ACS observations at shorter wavelengths. This 
wavelength baselin e of rest-frame UV to rest-fra me near-IR is optimal for constraining stellar age independent 



of dust extinction ( Gil de Paz fc Madorel 2002 ). We approach these multicolor data with several goals: (1) 



Understand the comparison of the stellar populations and color gradients of high-redshift AGN hosts with 
those of a non-active population of similar morphologies and mean redshift. Interpret this comparison in light 
of similar analyses at low redshift that find correlations between AGN strength and the presence of YSPs. 
(2) Characterize the morphologies of the hosts with diffraction-limited near-IR imaging and determine if the 
morphological classification is consistent with the stellar population distribution. (3) Compare the stellar 
populations of Type I and Type II AGN to shed light on the nature of AGN obscuration at high-redshift. 

To achieve these goals, we have chosen observational modes that deliver excellent spatial resolution 
(O'.'l ~ 600 — 800 pc) to separate bulge from disk and to distinguish nuclear star formation from the central 
AGN. Particularly for Type I AGN hosts, good resolution is necessary to subtract the point-like non-stellar 
continuum of the AGN itself, revealing the underlying host. At these redshifts and magnitudes (19 < R < 
23.5), stellar age determination via photometric SED fitting is preferable to spectroscopic analysis from the 
point of view of observing efficiency; AO resolution is necessary to study the bulge separately from the 
central point source and from the galactic disk, if either are present. 

We describe the data and analyze selection biases in the AGN sample in §2. In §3, we utilize multiple 
methods of AGN/starburst discrimination to demonstrate that the sample galaxies are indeed AGN. §4 
explains how annulus photometry and point source subtraction is performed on all imaging data. Section 
5 includes morphological analysis via multi-component GALFIT ( Peng et al. 2002h models of each band. 



with arguments for the presence or non-presence of a central point source in each band. In §6, we present 
B-i and V-K' color maps and multi-color surface brightness gradients for all AGN galaxies. Section 7 gives 
detailed stellar populations analysis of the color maps for the AGN sample. We make conclusions in §8 
and summarize in §9. We outline the method used to estimate the Keck Laser Guide Star AO (LGSAO) 
point spread function (PSF) in Appendix A. We give morphological and stellar populations analyses for a 
comparison sample of normal (non-active) galaxies in Appendix B and C. 



We assume a flat cosmology throughout, with Ha = 70 km/s/Mpc, ^Im = 0.3, and JIa = 0.7. 
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MULTIBAND DATA 



2.1. GOODS-South Optical, X-ray, and Ground-Based IR Data 



We investigate AGN in the Great Observatories Origins Deep Survey (GOODS) field South, or Chandra 
Deep Field South (CDF-S), an ideal field in the optical/IR both for depth and wavelength coverage. The 
wide wavelength baseline in this field enables the use of SED-fitting to estimate fundamental properties of 
the underlying stellar populations of galaxies. The deep, 1 Megasecond Chandra exposu re in this region has 



revealed hundreds of AGNs at 0.5 < z < 1.5 jPan et alJboOll . boOst iBarger et allboosh . HST ACS images 



of the GOODS fields (ICiavalisco et all 120041) in the F435W (B), F606W (V), F775W {i), and F850LP (z) 
filters are the most sensiti ve optical ima g es to date in these regions. We obtain spectroscop ic redshifts 



Giacconi et al 



(200 



in GOODS-S/CDF-S fr om ISzokolv et all (|2004l ) and photometric redshifts from COMBO-17 (|Wolf et al 
2004 : Zheng et al. 20041 ). In this paper, all source ID numbers are in the form "XID #" as reported in 



We also make use of deep Ks imaging in the CDFS using the ISAAC instrument mounted at the Antu 
Unit Telescope of the VLT. All sources are detected in this Ks imaging . The data reductio n and depths are 
described at the GOODS/ISAAC Final Data Release (V2.0) webpage |Retzlaff et al.lboOTi) . 



2.2. K' Band Adaptive Optics Imaging 



All K^-band adap t ive optics imagin g has been obtained as part of the CfAO Treasury Survey (CATS, 



Melbourne et al 



2005 



Koo et al.ll2007n . with three nights of laser and natural guide star imaging on the 



Keck II telescope. Table 1 lists pointings and tip/tilt star coordinates. These pointings were specifically 
chosen to optimize correction of the Chandra sources in the available field of view (FOV) of the adaptive 
optics system. Field centers were chosen to maximize the number of X-ray sources as well as serendipitous 
field galaxies in the FOV; the latter of which are used as galaxies for comparison in this study. The K images 
were obtained with the NIRC2 wide camera (40'.' x 40'.') with a pixel sampling of 0'.'04. This provides good 
sampling of the Keck telescope diffraction-limited PSF at 2.2 microns, whose full width at half maximum 
(FWHM) was 0'.'054. 

Imaging and reduction strategies are detailed in Melbourne et al. ( 20051 . 2007), with a cursory summary 
here. The tip/tilt star, varying in brightness between i? = 15 and i? = 17 for different fields, is included in 
the 40'.' X 40'.' corrected field of view to enable good temporal sampling of the PSF. The spatial variation of 
the PSF is captured with nightly observations of a globular or open cluster. For dim PSF stars {R < 16), 
we obtain 2 exposures of 2 coadds of 30 seconds at each dither position. The typical exposure of one hour 
consists of 30 dither positions marking random locations on a circle of radius 2'.' The reduction entails the 
creation of super-sky and fiat-field frames by median-combining the object frames and masking out sources. 
UKIRT IR standards were used to determine the photometric zeropoint. 



2.3. Sample Selection 

GOODS-South is an ideal field from which to choose targets for purely photometric stellar populations 
analysis, as deep ACS photometry is available in four bands, covering a wide wavelength baseline. The sample 
of AGN in this study is X-ray selected, using the catalog of sources extracted from the 1 Ms Chandra Deep 
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Field South (CDFS) observation as described in lGiacconi et alj (j2002l ). The flux limits of this final catalog 



for the soft (0.5 - 2.0 keV) and the hard (2-10 keV) bands are 5.5 x 10"" ergs s 



and 4.5 x 10 



-16 



ergs s ^ cm ^ 



, respectiv ely. Sources were extracted with a modified version of the SExtractor algorithm 



(jBertin fc Arnoutd 119961 ) with a S/N detection threshold of 2.1, where single-band detection was sufficient 
for inclusion in the c atalog. The extraction methods and accompanying selection biases are described in 
Iciacconi et al. ( 2002). Thro ughout this paper, th e hardness ratio is defined as HR = {H — S)/{H + R) as 
in Giacco ni et al.l (j2002l ) and lSzokolv et al.l (j200j), where H and S are the net instrument count rates in the 
hard (2-10 keV) and the soft (0.5-2 keV) bands, respectively. 

The Keck LGSAO system requires a point source in the field with i? < 17 to serve as a tip/tilt reference. 
Unfortunately, few stars are available in GOODS-South to use as tip/tilt references. Therefore, our sample 
of galaxies is limited by the number of available tip/tilt stars in the field. Fields with bright tip/tilt stars 
were observed first to maximize science return, as the quality of adaptive optics correction suffers for dim 
tip/tilt reference stars. Little selection bias is introduced by concentrating on targets with nearby bright 
stars. Of the 346 sources in the CDFS X-ray catalog, ~ 30 sources are present within the ACS fields, possess 
suitable tip /tilt sta rs nearby, and have photometric or spectroscopic redshifts measured (jZheng et al.ll2004 : 
Szokolv et al1l2004h . 



We also impose an R band magnitude cutoff for each galaxy of 23.5. X-ray sources with optical coun- 
terparts dimmer than this value do not possess sufficient signal to perform stellar populations analysis in 
the outer regions of the host (r > I'.'O) in any band. Of the 12 AGN sources imaged thus far by the CATS 
survey, this magnitude cut removes three sources: XID 629 (HR = -1.00, R = 25.26), XID 532 (HR = -0.06, 
R = 24.59), and XID 61 (HR = -0.45, R = 23.98). XID 629 appears in the ACS bands as a compact, 
asymmetric structure on the order of 1-2 kpc in half-light radius (assuming z > 0.5), with a soft X-ray 
spectrum that would classify it as a starbursting galaxy or Type I AGN. XID 629 does not appear in the 
LGSAO K' image. XID 532 appears in the ACS bands as a powerful point source surrounded by a streak of 
compact, asymmetric star formation and would be classified as a Type II object or starbursting galaxy from 
hardness ratio alone. XID 61 is a pure optical point source, which when taken with the hardness ratio of 
-0.45, classifies it as a likely Ty pe I QSO. XID 61 is a dim source located 4'.' away from a bright star, which 
prevented good spectroscopy in ISzokolv et al.l (|2004r ). although a photometric redshift is available. 



Apart from the R = 23.5 magnitude cut, we also remove sources that are clearly dominated by a point 
source. It is challenging to isolate and analyze the hosts of these compact or point-source dominated sources. 
XID 30 (HR = -0.53, R = 20.42, z = 0.837), which has been cla ssified as a Type I b road-line AGN (BLAGN) 
from optical spectra and as a Type I QSO from X-ray spectra (jSzokolv et al.ll2004l ). appears as a pure point 
source in the ACS bands. A detailed two-component fit to the 2D image (see §5.1 for more information 
on GALFIT) reveals that the central point source is more than 3 times brighter in i-band than the total 
luminosity of any underlying host. 



Properties of the remaining sample are enumerated in §3.1. 



3. AGN AND COMPARISON SAMPLES 



Following the magnitude and morphology cuts detailed in §2.3, eight AGN remain in the sample. Co- 
ordinates, redshifts, and broad photometric properties are given for these sources in Table 2. 
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3.1. Verifying AGN identification 



Identifying AGN in high-redshift surveys of any wavelength has been historicaUy difficult due to their 
widely variant morphologies and spectral properties. In particular, galaxies that are undergoing vigorous star 
formation may display X-ray or IR luminosities that mimic galaxies that harbor nuclear activity. Multiple 
historical methods of distinguishing AGN from starbursts exist, including (1) total X-ray luminosity or flux 
ratio of x-ray to optical; (2) particular optical narrow-line ratios like [OIII]X5007/Hf3 and [N 1 1] X6583/ Ha; 
(3) flux ratios in the near- and mid-IR; and (4) presence of broad-lines in the optical spectrum. In this 
section, we use combinations of evidence from all four categories above to demonstrate that all the objects 
are truly AGN. These methods of discrimination are described in more detail in the following paragraphs. 
Shown in Table 3 are X-ray, infrared, and optical line luminosities and ratios for each of the sample galaxies. 



Sources with total X-ray luminosity Lx. 



0.5-lOfcey 



greater than 10"^ ^ ergs s ^ are likely to be AGN, as 



few l ocal starbursting galaxies have luminosities greater than this value (jRosati et al.ll2002t [Alexander et al 
20011 ). T he X-ray c lassification shown in Table 2 is derived from a two-dimensional scheme originally pre- 



sented in lSzokoly et al. (.2004. ) . in which hardness ratio purely determined AGN Type and X-ray luminosity 
purely determines strength. Objects with hardness ratio less than -0.2 are Type I objects; those with 
hardness ratios greater than -0.2 are Type II objects. X-ray sources are classified as "AGN" if Iff*^ ergs 
s"^ < Lx,o.5-iokeV < 10'*'^ ergs s"^ for Type I objects and lO''^ ergs s"^ < ix,o.5-iofcey < 10'*'' ergs s~* 
for Type II objects, or clas sified as "QSOs" if Lx,o.5-iokeV > ergs s^^. Local hard X-ray surveys find 
very good correspondence (jTueller et al.ll2008[ ) between X-ray typing by hardness ratio or neutr al hydrogen 
column density and traditional typing by characterizing widths of lines in the optical spectra as in lOsterbrock 
(Il989h . 



The X-ray-to-optical flu x ratio is a good d i scrim inator between AGN and normal galaxies. The original 



formulation for local AGN in Maccaca: 



Maccacai 


ro et al. (1988) 


difled by 


Szokolv et al. 



( 1988 ) is a f unction of the 0.3-3.0 keV soft band flux and the 



V-band magnitude. Since modifled bv lSzokoly et al.l (|200J) for higher redshift objects, the X-ray-to-optical 
flux ratio (X/0 hereafter) in Table 3 uses the R-band magnitude and the 0.5-2.0 keV Chandra soft band 
flux. AGN at higher redshift tend to satisfy —1.4 < X/O < 1.1 and normal galaxies have —2.2 < X/O < 0.0. 
Table 3 plots values of the X/O flux ratio for all systems, including a correction in the Type II systems for 
soft-band absorption by dust. 

An add it ional tool we utilize is the 24 micron to 8 micron flux ratio (24/8 hereafter), as discussed in 



Brand et al.l (|2006f ). At high redshift a histogram of 24/8 ratios for a mid-IR selected sample is bimodal, 
centered about 24/8 = and 0.5. Sources clustered about are identified with AGN that produce mid-IR 
SEDs similar to power laws; those clustered about 24/8 = 0.5 are identified with starbursting galaxies that 
possess SEDs characteristic of warm dust. 



3.1.1. Type II AGN 



XID 155, 56 and 266 are identified as Type II AGN by hardness ratio and X-ray luminosity. These 
systems also have X/O ratios between -0.1 and -0.7, which places them in the AGN range. XID 56 
displays optical spe ctral signatures of A GN, including [Ne V]. [OIII] A5007 and Hf] line intensities are 



taken directly from iNetzer et al 



( 2006h . with an upper limit on iJ/3 estimated for XID266. The [OIII] 
X5OO7/H0 line ratio is between 4 an d 8 for all systems, placing them in the upper half of a BPT diagram 
([Baldwin. Phillips, fc Terlevichlll98ll ). Although the X-ray properties of these systems clearly classify them 
as AGN, their 24/8 flux ratios are all greater than 0.3, which implies that their mid-IR SED includes some 
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contribution from a starburst. 



3.1.2. Type I BLAGN 



XID 15 and 32 are broad-line AGN based on their optical spectra. Their X-ray luminosities {Lx,o.5-iokeV = 
63 X 10''^ ergs s~^ and 11.1 x lO'*^ ergs s~^, respectively), X/0 ratios (0.2 and -0.1), and hardness ratios 
additionally classi fy them as X- r ay Ty pe I AGN. We estimate the [OIII] A5007 line luminosity of XID 32 
with spectra from lSzokolv et aL (2004). XID 15, however, is at too high of a rcdshift (z = 1.23) to measure 
the [OIII] A5007 strength with an o ptical spectrum, so we estimate its strength via a log-linear calibration to 
hard x-ray luminosity presented in 



Netzer et all (|2006l ) for Type I AGN and QSOs. It appears from Figure 



2 of that work that the precision of this calibration is half a decade. 



3.1.3. Type I AGN with no optical spectral signatures 



XID 536, 594, and 83 are X-ray sources that have no optical signatures of the presence of an AGN. 
XID 594 is classified as a Type I AGN by its X-ray luminosity of 2.44 x 10"^^ ergs s~^, its X/0 ratio of -1.1, 
its 24/8 ratio of -0.2, and its hardness ratio of HR = -1.0. XID 83 is similarly classified as a Type I AGN 
by these same properties {Lx.o.h-iokeV = 4.4 x 10^^ ergs s~\ X/ = -0.1, and 24/8 of -0.18). XID 83's 
redshift is derived from a photometric analysis (jZheng et al.ll2004^ with a 2a confidence interval of 1.68 to 
1.95. A true redshift of less than 0.5 would imply that the true 0.5-10 keV X-ray luminosity were less than 



10"*^ ergs s ^, but such a scenario is unlikely for this red {Rab 



K' 



= 5.81) source. XID 536 has an 
X-ray luminosity and X/0 ratio below the AGN threshold {Lx,o.5-iokeV — 3.0 x 10^^ ergs s'^ and X/0 
= -2.3) and is thus classified as a galaxy in Table 2. This X-ray luminosity, however, would still require a 
substantial starburst to explain. As we show in §8, the stellar populations of this object are not consistent 
with a starburst. In addition, the 24/8 flux ratio of -0.40 is quite low for a starburst, but can be better 
accounted for with low AGN activity. We surmise that XID 536 is a weak AGN and therefore include it in 
the sample as such. Further, its neutral hydrogen column density of less than 3.0 x 10^^ cm~^ argues for an 
unobscured. Type I source. 



We estimate the [OIII] A5007 line luminosity of XID 594 with spectra from lSzokolv et aL ( 2004) and 



show the results in Table 3. For XID 536 and 83, sources that do not have optical spectra in lSzokolv et al 



(2004), we estimate the [OIII] A5007 line luminosity using the X-ray luminosity calibration of 



Netzer et al 



20061 ) as applied on XID 32. 



3.2. The Comparison Sample 

Stellar populations analysis has been shown to be prone to systematic errors in various cases. It is 
useful to compare the stellar populations statistics of the sample of interest to a control sample at the same 
redshift. In this way, the systematic errors present in the stellar populations fitting will be present in both 
samples and the relative differences may be investigated. 

Nine non-active galaxies have been chosen from the 40'.' x 40'.' AO corrected field of view to serve as 
a comparison sample. All selection criteria imposed on the AGN sample have also been imposed on this 
sample. In addition, the galaxies were selected to have a similar distribution of integrated K' apparent 
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magnitudes. This necessitated choosing the brightest galaxies in the adaptive optics K' images to match 
the AGN. The system coordinates, R-band magnitudes, and redshifts are given in Table 2. The mean R 
and K' (AB) magnitudes for the comparison sample are 22.5 and 20.1, respectively, as compared to 21.65 
and 19.3 for the AGN sample. The mean R-K' color for the comparison galaxies is 1.40, as compared to 
2.35 for the AGN sample. The mean redshift is 0.724 for the comparison sample and 0.836 for the AGN 
sample. It is expected that the mean redshifts should be similar, and the small difference seen here does 
not affect the broad conclusions of the paper. K' is a rough proxy for stellar mass, and AGN have been 
shown to preferentially lie in massive galaxies (e.g., Kauffmann 2003). Thus, choosing the most K'-bright 
(or most massive) non-active systems to compare to an active sample is suitable. The mean K' magnitude 
of the AGN sample is still 0.8 magnitudes brighter than that of the comparison sample, which reflects our 
inability to find sufficiently massive galaxies for the comparison sample in the limited survey field of view. 
These minor differences will be taken into account in later analysis. 



4. PHOTOMETRY 

4.1. Estimation of LGSAO Point Spread Function 

Adaptive Optics systems tend to produce point spread functions (PSFs) that are the sum of two com- 
ponents, a diffraction-limited core and a wide halo with FWHM roughly the size of the seeing disk. The 
Strehl is a measure of the amount of power in the core of the PSF relative to a perfect PSF produced by 
the telescope and instrument. The PSF in AO systems is also spatially and temporally variant, complicat- 
ing precise photometry on spatial scales comparable to the diffraction limit. The PSF must be well-known 
to permit the photometric accuracies we seek for the stellar populations analysis in this work. We use a 
semiempirical method of determining the PSF, which utilizes an LGSAO image of a globul ar or open cluster 



taken once throughout night and the PSF of the tip/tilt star within the corrected FOV (jSteinbring et al 



20041 l2005h . We present this method and the nuances of its application to this work in Appendix A 



4.2. Subtraction of Central Point Source 

AGN and QSOs tend to have contamination by non-thermal continua in the central regions, due to the 
accretion disk. This feature is expected to be unresolved and bluish. Unfortunately, the feature we search for 
(circumnuclear starbursting) is typically less than 1 kpc in extent and is likely hidden within the non-thermal 
contamination. The Keck diffraction limit at K is 54 mas, which corresponds to ~ 400 pc at a redshift of 
0.8. Thus, the principal strategy is to estimate stellar populations statistics surrounding the central point 
source (r > 1 kpc) and ignore the stellar populations of the central point source itself. 

However, central unresolved point sources contribute scattered light to the outer regions of interest, so 
it is necessary to subtract the central point source before attempting aperture photometry. We determine 
whether or not a point source is present in a particular band with GALFIT. We first fit the most radially 
extended components of the galaxy with two sersic profiles (bulge and disk) on the i-band image. Then, we 
introduce a Gaussian component with a range of FWHM between — 3 pixels, or less than 100 mas, and 
refit all other parameters (including the Gaussian magnitude) . If the variation in between the zero- width 
Gaussian case and the 3 pixel Gaussian case is not significantly different at the 5-sigma level according to 
an F-test, the central peak is not judged to be a point source. If the central peak is judged to contain a 
point-like component, it is subtracted from the images before annulus photometry is performed. 
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By the criterium above, full subtraction of a central point source was necessary for two sources, XID 
15 and XID 32, both Type I sources and BLAGN by optical spectrum. To indicate that the results of the 
point source subtraction are sensible, we plot the best fit SED of the removed poin t sourc e for each source 



in Figure 1. These SEDs are plotted against quasar color models from lFrancis et al.l (j200Cll ). It appears that 
a power-law SED {F^, oc with a small infrared emission component is a sensible fit to the point source 

SED of XID 32. The RMS error to this fit is 0.06 mags, which is on the order of the expected point source 
subtraction error. A power-law SED with a small dust absorption component, E{B — V) = 0.3, fits the SED 
of XID 15 to an RMS error of 0.07 mags, which is also on the order of the expected error. 

Subtraction of central point sources is necessary to analyze host stellar populations in AGN and QSOs, 
but the point source SED in itself gives information about the extinction seen by the nuclear emission. For 
those systems that lack a central point source in a particular band according to the criterion above, we derive 
an upper limit on the magnitude of a point source by subtracting the brightest possible PSF from the core 
of the galaxy that does not invert the radial profile in the central regions. Typical galaxies are not expected 
to have "brightness holes" in the nucleus, so the brightest PSF that subtracts to produce a flat radial profile 
is an upper limit on the magnitude of a central point source. 



4.3. Aperture Corrections for Annulus Photometry 

Knowledge of the intrinsic brightness profile of a galaxy is limited by PSF broadening in this study. 
Here, we compose sersic-type models of the intrinsic distribution and convolve them with the known PSF in 
an attempt to match the data. After finding the best model, we compute the annular "aperture corrections" 
that account for PSF broadening. These corrections are then applied to the photometry measured on the 
reduced (but unprocessed) data. We use GALFIT to perform this model-fitting, dividing the unconvolved 
model by the convolved model to yield the aperture corrections. 

We perform annulus photometry over 10 successive, non-overlapping circular annuli with increasing 
radial widths of 0'.'15. This gives a central core aperture with a diameter of 300 mas. 

This method is illustrated in Figure 2 as a series of radial profiles taken at intermediate steps during the 
calculation of the aperture corrections for XID 56 (K'). Shown are the raw radial profile, the best GALFIT 
estimate of the intrinsic light distribution, and the corrected annular photometry. 



4.4. Errors in Annulus Photometry 

We choose a wide central aperture because the variance of the 300 mas encircled energy across the 
diameter is typically 2-3 times lower than the Strehl variance across the field, especially for our high airmass 
fields. The errors in our PSF estimation method are manifest on small spatial scales, as they are due to 
misestimates of the isokinetic angle, so only the first radial point is affected in the annulus photometry. 
We do not attempt to draw conclusions about morphology on spatial scales on the order of the Keck K 
diffraction limit (^^ 54 mas). 

Errors are calculated in each annulus by summing in quadrature the following error sources: (1) photon 
noise, (2) read noise, (3) dark current noise, (4) sky noise, (5) PSF subtraction error, and (6) aperture 
correction error due to uncertainty of the PSF in K. Errors 2-4 are computed by the root mean square of 
an empty sky region. Source (5) refers to the errors in estimating the photometry of the host galaxy caused 
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subtracting a central point spread function. In point source-dominated galaxies (like Type I QSOs), this 
error increases beyond 0.3 magnitudes in the central regions. It is estimated with a Monte Carlo technique 
of simulating the accuracy of PSF recovery for a particular galaxy. A central core PSF is added to the host, 
which is first obtained by subtracting the original central PSF. GALFIT is then used to refit the galaxy 
components plus PSF, and the best-fit PSF magnitude recovered. The GALFIT initial conditions are varied 
and the 1-sigma standard deviation of recovered PSF magnitudes about the true simulated magnitude is 
calculated. This error source is strongly centered in the first annular core and drops sharply afterward. 

Error source (6) is computed with a separate Monte Carlo technique to determine the effect of PSF 
misestimation on final host galaxy photometry. As discussed above, the most uncertain factor in the PSF 
estimation is the isokinetic angle, which strongly affects the encircled energy in the first few annular rings. 
This error plays into the host galaxy photometry through the aperture corrections. We compute the error 
on the aperture corrections by simulating 100 PSFs with a distribution of isokinetic kernels convolved (i.e., 
Gaussian kernels with varying FWHM). The aperture corrections are computed in each case and the one- 
sigma errors are derived from the resulting distributions for each annular ring. 



5. MORPHOLOGICAL PROPERTIES OF INDIVIDUAL SOURCES 

We now describe the morphological properties of each of the sources in the AGN and comparison 
samples. This section attempts to answer the following questions for each AGN host: 1. What are the broad 
morphological properties of the system? 2. What limits can be placed on the presence of a point source in 
each band? The answer to question (2) may influence our interpretation of each system as an obscured or 
unobscured AGN. Specifically, the SED of such a point source can be used to constrain the extinction seen 
by the central emission. In this section, we answer these questions for each AGN host. Figures 3-10 show 
B-i and V-K' colormaps for each AGN source to explore trends in color gradients, as well as tiled imagery 
and model fits. A summary is presented at the end of the section. We display image tiles, GALFIT fits, 
model residuals, and morphological analyses for the normal, comparison sample in Appendix B. 

Colormaps are generated with the following process. For a particular galaxy, images from each of two 
bands are convolved with the PSF of the other band to match resolutions. Then, to preserve signal in 
the noisy wings of the galaxy, the images are both convolved with a Gaussian kernel with spatially variant 
FWHM. The FWHM increases linearly in a radial pattern, from zero in the center to I'.'O at a radius of 
0'.'75. Thus, good spatial resolution in the colormaps is only preserved in the central square arcsecond. The 
V-K' colormaps utilize near-IR data from LGSAO imaging in the central 2'.'0 and deep ISAAC imaging in 
the surrounding regions. 

5.1. Type II AGN 
5.1.1. XID 56 



XID56 is a clear merging system, with central radially symmetric emission closely surrounded by blue 
spiral arms. Several point sources are seen in the central 10 kpc, including one blue source offset from 
the center that aligns with a line of star formation in the northeast. This blue lin e of star formation ma y 



alternatively be the optical counterpart of a jet associated with this point source (jMelbourne et al 



200. 



i). 



The GALFIT fit to the K' band image possesses a large, extended spherical structure lying beneath the 
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many features in this merger. Both the B-i and V-K' colormaps reveal this component to be redder than 
the bluer spiral arms and knot-like structures. The central structure is not well-fit by a point source in any 
band. 



5.1.2. XID 155 

XID 155 has two major components: A large vertical disk-like structure lying beneath a compact, central 
source is visible in the V-K' colormap shown in Figure 4. The central structure is not well-fit by a point 
source in any band. The underlying inclined disk structure, seen in both B-i and V-K' colormaps, is bluer 
than the central source. This appears to be a disk with a significant bulge component. 



5.1.3. XID 266 

XID 266 is an inclined star forming disk with an appreciable central bulge component. Small tidal tails 
and star-forming knots on the North and South ends of the disk, as well as a possible blue companion 1" to 
the NW, indicate that minor interactions are probable. The disk star formation is manifest in multiple star 
forming knots that arc pinned to two major spiral arms. Judging from the connection between the spiral 
arms and the bulge, there appears to be a central bar as well. Both B-i and V-K' colormaps indicate strong 
negative color gradients, or redder in the bulge and bluer in the outer regions, which is consistent with disks 
with large bulges. Central point source models do not fit well to the bulge in any bands. 



5.2. Type I AGN 

5.2.1. XID 83 

XID 83 is a spherically symmetric object with some hints of a disk and an appreciable spheroidal 
component. Small spiral arms are apparent in the log-scale images and are seen clearly in model-subtracted 
residuals in all bands. The nearby blue structure in the SE corner was masked for measuring photometry. 
The B-i colormap indicates that the central core is bluer than the underlying spheroidal component. The 
V-K' colormap indicates this positive color gradient as well, but the overall system is quite red in V-K' and 
is very bright in K' out to a radius of 1.5 — 2'.'0. Red observed frame V-K' and bluer B-i colors are expected 
for such a high-redshift object (z = 1.76). The very red, dominating spheroidal component underlying a 
bluer core may be indicative of a large elliptical with central star formation. Central point sources models 
do not fit well to the central core in any optical band, although a small [R,, < 1 kpc) knot is possibly present 
in K'-band. 



5.2.2. XID 536 



XID 536 is a clear merger of two spheroidal galaxy cores, previously presented in lMelbourne et al.l (|2005f ). 
Neither central core is a point source in any band, although individual sersic profiles fit well. The underlying 
spheroid is well-fit by a single, radially symmetric sersic profile with no disk-like signatures. Present in 
the optical model-subtracted residuals are tidal features and shells, which indicates kinematic interaction 
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between the two cores. These features may also be remnants of spiral arms, although they do not display 
knot-like star formation. The B-i colormap indicates that the underlying spheroid is redder than the Western 
(right) core and bluer than the Eastern (left) core. Both cores are bluer in V-K' than the surrounding ~ 1", 
but about the same color as the outer regions (1" < r < 2"). 



5.2.3. XID 594 

XID 594 is a spheroidal galaxy with a bright central core and no disk- like features. All bands display 
smooth morphology in the model-subtracted residuals. Both the B-i and V-K' colormaps indicate that the 
central core is very red, while the underlying spheroid is slightly bluer, but still red {B—i ~ 4.0, V—K' ~ 3.5). 
The central core is larger than a point source in all bands. 



5.2.1 XID 15 

XID 15 is a broad-line AGN (BLAGN) with a clear point source present at the center in all bands. 
Surrounding the central, blue point source is an inclined ring of blue star formation (axis ratio = 0.77). 
There do not appear to be any other significant morphological features. The SED of the centra l point source 



i s con sistent with a dust-reddened QSO {E{B — V) ^ 0.25 by comparison with Figure 1 of iFrancis et al 
(j2000l )) and a spectral index a — —0.3 (i.e.. Si, oc z^"). 



5.2.5. XID 32 



XID 32 is a broad-line AGN with a clear, blue point source in all bands. The SED of the point source is 
consistent with an unreddened QSO with some infrared emission (e.g., see Figure 2 of iFrancis et~aLl |2000[) ) 
and a spectral index a = —0.3. Unlike XID 15, this BLAGN sits on a spheroidal component which is redder 
than the point source. This component is eight times dimmer than the point source in observed B band and 
three times brighter in observed z-band. 



5.3. Summary of Morphological Properties 

The AGN in this sample represent a range of morphologies, including irregulars/mergers, normal star 
forming disks with bulges, and ellipticals. All of the AGN sources except one (XID 15) have an appreciable 
bulge or spheroidal component. Although there appears to be no preferred morphology for AGN hosts, a 
correlation between morphology and AGN Type is supported. All of the Type II AGN either possess disks 
(XID 155 and XID 266) or are irregular mergers (XID 56). One of the broad-line Type I AGN is a star 
forming ring, which indicates disk-like morphology, but four of five of the Type I AGN are ellipticals (or 
merging ellipticals, in the case of XID 536). This is consistent with a previous study of morphology of 
AGN hosts at high rcdshift in the Extended Groth Strip (EGS) survey bv lPierce et all (120071 ). The authors 
utilized Gini-M20 non-parametric analysis on Chandra sources in EGS, finding that the median hardness 
ratio of E/SO/Sa galaxies is -0.46 (Type I unobscured) and the median hardness ratio of Sc/Irr/d hosts is 
0.55 (Type II obscured). 

The normal galaxies also represent a range of morphologies, including 6/9 disks, spirals and irregulars 
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and 3/9 ellipticals. Ellipticals are 50% more common in the AGN sample (4/8) than in the normal sample 
(3/9), but this is not significant for this small sample. No conclusions about the comparison between AGN 
hosts and nonactivc samples can be drawn from the morphologies alone. However, the sample of normal 
galaxies will be helpful in the following stellar populations analysis in that they may reveal systematic errors 
in this fitting process. 



6. STELLAR POPULATIONS ANALYSIS 

For eac h annulus, aperture-correc ted photometry is available in observed B, V, i, z, and K. These are 



compared to lBruzual fc Charloti (j2003f ) tau models at the redshift of interest. Bruzual & Chariot tau models 
predict the integrated spectrum of a stellar population at multiple timepoint s thou gh its existence, with 
a star formation rate that depends exponentially on time. We use a ISaloeteij (|l955l) initial mass function 



(IMF). These models are obtained at high redshift by blueshifting the filter pass functions of the requested 
filters as well as blueshifting the input AO spectrum (as in Melbourne et al 2008). All model grids are 
computed over a range of three variables: (1) Age since star formation begins; (2) tau parameter (e-folding 
time) in years; and (3) the extinction resulting from the dust mix ed with the stars. A GN have been shown 
to have solar or super-solar metallicities at a range of redshifts ( Groves et al. 2006[ ). so we assume solar 



metallicity. The Bruzual & Chariot models are interpolated for 210 ages, 16 values of dust extinction, and 
7 values of the tau age. The best model is chosen in a least-squares sense, using the measured photometry 
with errors as derived above. Since the number of interpolate models is discrete, and uncomfortably low in 
the case of the tau parameter, a cubic interpolation is used on the final model values to obtain finer accuracy. 

The errors on the stellar population parameters are computed with a further Monte Carlo simulation 
utilizing knowledge of the photometry errors. For each ring and source, the least-squares fitting code is 
repeated on 100 manifestations of the photometry, varied according to the known errors. The resulting 
distributions of age, tau, and dust are recorded and used to derive the one-sigma errors for each parameter. 

The fractions of the population (by mass) are computed from the best-fit tau model for three different 
age ranges: Young (age < 100 Myr), Intermediate (100 Myr < age < 1 Gyr), and Old (age > 1 Gyr). All 
references to "young," "intermediate," and "old" refer to these age ranges. These fractions are calculated 
with integrals over an exponentially declining star formation rate. These fractions by mass do not refer to 
the fraction of a population after aging and undergoing total mass loss from the death of young stars; they 
refer to the fraction of all stars formed by mass in a particular time frame since star formation began. 

We stress that since all models are searched for every reasonable value of age, tau parameter, and dust 
extinction, the error bars we generate with Monte Carlo simulations are a good indication of the degeneracies 
in the model for a particular SED. Some SEDs may be fit by vastly different stellar populations models to 
equivalent in the presence of equivalent photometric error, while other SEDs are much less degenerate. For 
example, SEDs with extremely red (optical - K') colors may rule out all dust-less old or young populations, 
as these have a limit to their redness; in this case, the well-known age/dust degeneracy is mitigated. As such, 
one can determine if model degeneracy is playing a role by observing the size of the error bars relative to 
the fraction in each population; if the error bars in each population overlap, little conclusion can be drawn. 
We also point out that systematic errors in the models and insufficiencies in the type of model chosen (tau 
models with fixed solar metallicity) are present for both the AGN and comparison sample. The principal 
conclusions of this work are drawn from relative differences between subgroups of this sample. 



To illustrate the degeneracies of the models, Figure 11 plots photometry and best-fit spectra for two 
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cases: (1) the SED from the fourth ring of XID 155, which is representative of non-degenerate cases and 
(2) the SED from the fifth ring of XID266, which is representative of more degenerate cases. Each plotted 
spectrum represents a fit to the photometry when stellar age has been fixed, but dust extinction and the tau 
parameter are minimized. The solid lines in both cases represent the best-fit models, fixing no parameters 
(except metallicity) . Notice that the particular SED of the degenerate case may fit ages between 500 Myr 
and 4 Gyr, but the SED in the non-degenerate case does not fit the 500 Myr case. In both SEDs, young 
stellar ages (< 100 Myr) are ruled out by the near-IR point. Note that in Figures 13 and 14, the error bars 
on the best-fit model parameters refiect the level of model degeneracy apparent in Figure 11. 

The stellar populations fractions as described above are plotted for each host galaxy as a function of 

radius for the AGN sample in Figures 12-19. Also shown in each figure is a plot of radial profiles for each 
band. Appendix C contains the stellar population fraction plots, radial profiles, and short discussions for 
the normal, comparison sample galaxies. 

6.1. Type II AGN 

6.1.1. XID 56 

The merging system XID 56 displays predominantly young stellar populations in the core, older stellar 
populations in the range 0'.'4 < r < O'.'T, and intermediate age populations in the outer regions. The 
outer photometry is likely from the linear feature in the Northeast, so the intermediate ages refer to this 
structure. The underlying red spheroidal structure discussed in the previous section is composed of older 
stellar populations. 

6.1.2. XID 155 

The inclined disk XID 155 possesses a central core composed of dusty younger stellar populations. The 
outer disk (r > 0'.'2) is largely old, with 4-6% of an intermediate age population by mass and negligible young 
populations. 

6.1.3. XID 266 

The large star- forming disk XID 266 displays a prominent negative color gradient, with a red core and 
blue outer regions. The central bulge appears to be a mixture of all ages, but is dominated by older stellar 
populations. The star forming disk (r > 0'.'2) is a mixture of intermediate aged populations. The fraction of 
older population associated with the spheroidal component decreases with radius. 

6.2. Type I AGN 

6.2.1. XID 83 

The red spheroid XID 83 with small central spiral arms displays intermediate aged populations in the 
central 0'.'6 and a mixture of old and intermediate aged populations beyond this radius. This suggests that 
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the central spiral arms are the remnant of a recent star formation event. 

6.2.2. XID 536 

The elliptical merger XID 536 is dominated by old stellar populations at all radii, although intermediate 
aged populations make a strong contribution between 0'/2 and 0'/7. Both the two galaxy cores and the streaks 
that appear in the optical residuals fall in this region and may represent these populations. 

6.2.3. XID 594 

The smooth elliptical host XID 594 is dominated by old stellar populations at all radii greater than 0'.'15. 
The central core is mostly composed of intermediate aged populations, although the error bars indicate that 
the models may be degenerate here. 

6.2.4. XID 15 

The broad-line AGN XID 15 is largely composed of intermediate-aged stellar populations at all radii. 
The radial profiles shown in the left panel of Figure 18 and used to reconstruct the stellar populations 
statistics plotted in the right panel have been measured following a subtraction of the central point source 
in all bands, which is the reason that the colors of the central core do not resemble the colors elsewhere at 
higher radii. The point source contribution is estimated with a simultaneous sersic/PSF fit using GALFIT. 
As a result, the residuals are forced to resemble a galaxy-like sersic profile. There is a peak in the fraction of 
the youngest stellar population at a radius of ~ 0'.'7, which corresponds to the location of the star forming 
ring. 

6.2.5. XID 32 

The broad-line AGN host XID 32 is dominated by older stellar populations at all radii. Similar to 
XID 15, the radial profiles and reconstructed stellar populations statistics have been measured following a 
subtraction of the central point source in all bands using GALFIT. The old stellar populations arc consistent 
with the redness of the spheroidal component underlying the central point source. The contribution from 
young stellar populations is less than 0.3% for r > Cf.'2. 

6.3. Summary of Stellar Populations Analysis 

A summary of the dominant stellar populations in the inner (r < 0'.'3) and outer regions (r > 0'.'7) is 
shown for all galax;ies in Table 4. A number of observations about the stellar populations of the normal 
sample inform the comparison to the AGN sample. First, young and intermediate age populations are rarely 
seen as dominant in the outer regions of the systems (r > 3 — 6 kpc, depending on redshift). Those galaxies 

with this property tend to be spirals, disks, and irregulars (Norm3 and Norm6). Second, young populations 
are only seen as dominant in the central regions of spirals, disks, and irregulars (Norm2 and Norm5). Third, 



- 16 - 



ellipticals are more likely to have old stellar populations at all radii compared to spirals, disks, and irregulars 
(Norml, Norm4, and Norm?). Finally, as these results are expected for normal galaxies, they lend credence 
to the stellar populations fitting method. 

The stellar populations of the AGN hosts are consistent with their morphological classification. All 
three observations made in the above paragraph can be made separately for the AGN population, without 
considering AGN Type (but provided that central point sources resulting from non-stellar continua are 
subtracted off). In fact, for every AGN host it is possible to find a rough analogue in the normal sample 
with respect to broad morphological and stellar population characteristics. This being stated, the correlations 
seen in §5 between morphology and AGN Type are manifest as correlations between the presence of young 
and intermediate stellar populations and AGN Type. A majority of Type I AGN (4/5) have old stellar 
populations in the outer regions. All of these 4 galaxies are ellipticals. A majority of Type II AGN (2/3) 
have intermediate stellar populations in the outer regions, although this is not significant due to small number 
statistics. 

The strongest summarizing observation that can be drawn from this small sample is that the Type II 
AGN representatives are more likely to be spirals, disks, or irregulars and the Type I AGN representatives 
are more likely to be ellipticals according to corroborating evidence from both morphology and stellar 
populations. 



CONCLUSIONS AND DISCUSSION 



7.1. Stellar Populations of AGN in Local Universe 



The stellar populatio ns of this z ^ 1 populati on of AGN may be compared to those at lower redshift. 
A study of SDSS AGN in iKauffmann et al.l (j2003l ). selected via optical spectroscopic line ratios find that 
narrow-line AGN (Type II, 0.02 < z < 0.3) of all luminosities reside almost exclusively in massive galaxies 
and have distributions of sizes, stellar surface mass densities and concentrations that are similar to those 
of ordinary early-type galaxies in our sample. Locally, the host galaxies of low-activity AGN have stellar 
populations similar to normal early types, AGN activity/strength being measured with the line luminosity 
of the [QUI] line0 The hosts of high-activity AGN have much younger mean stellar ages. The young stars 
arc not pref erentially located near th e nucleus of the galaxy, but are spread out over scales of at least several 
kiloparsecs. IKauffmann et al.l (j2003l ) also examines the stellar populations of the host galaxies of a sample 
of broad-line AGN and conclude that there is no significant difference in stellar content between Type 2 
SevfertI (|l943h hosts and quasars (QSOs) with the same [OIII] luminosity and redshift. This establishes that 
a young stellar population is a general property of AGN with high [OIII] luminosities. Strong AGN have 
stellar populations similar to late type galaxies (young ages). 

Although the correlation bet ween AGN strength measured with L^qiii] ^-i^d f^ic hard X-ray luminosity 
is good at low and high redshift ( Netzer et al. 2006[ ). the correlation between AG N Type determined with 
optical spectroscopic methods and X-ray hardness ratio is less clear. The AGN in IKauffmann et all |2003h 
are predominantly Type II AG N, determined by narrow-line ratios in optical spectra. The traditional typing 
method used for decades (e.g.. IOsterbrocklll989l ) requires the presence of optical broad lines to flag a Type 



'^AGN strength is independent of AGN Type. Strong AGN are defined as log i[o//7] > 7 in units of bolometric solar 
luminosities (Lq = 4 x 10^^ ergs s"-*^). 
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I AGN. The X-ray Type classification relies primarily on the X-ray hardness ratio (HR) and the neutral 
hydrogen column density that can be computed from the HR through models. This typing confusion has 
been addressed with deep hard X-ray surveys, like the Swift BAT survey (jTueller et al. 2008.1 . This survey of 
^ 100 local AGN finds that the correspondence between hard X-ray typing and optical typing is good (94% 
correspondence) when archival optical spectra of X-ray AGN are reexamined for broad-lines. We continue 
with the knowledge that X-ray AGN Types and strengths are comparable to Types and strengths derived 
primarily from optical spectra. 



7.2. Comparison to Local Results 



As lKaufFmann et al.l (|2003l ) divides the local Type II sample into two categories (strong and weak), we 
divide the high redshift sample into four categories (including Type I AGN): Strong Type I AGN, Weak 
Type I AGN, Strong Type II AGN, and Weak Type II AGN. The distinction between Type I and Type II is 
made at a neutral hydrogen c olumn density of 1 ^^ cni ~^. The division between weak AGN and strong AGN 
is made at log L^oiii] = 7.5. Kauffmann et al. ( 2003h uses a division of log L^qiii] = 7.0, but we require 



a higher division to include 3 galaxies in the weak sample rather than 2 (XID 32 is in the Strong Type I 
sample if the sample is divided here). Our sample of 8 AGN hosts with 0.4 < z < 1.76 has morphologies. 



stellar population distributions, and color distributions that are consistent with findings of iKaufFmann et al 



(|2003l) at lower redshift. 



Type I AGN in this sample are not dominated by young (age less than 100 Myr) populations at any 
radii. Type I AGN may be dominated by intermediate aged stellar populations, but these are restricted to 
the central (r < 3 — 6 kpc) regions. Overall, 5 of 5 strong AGN (both Type I and II) hosts are dominated by 
young or intermediate aged populations at some radius, while 1 of 3 weak AGN (both Type I) are dominated 
by young or intermediate AGN at some radius. To enable comparison across the AGN Type dimension as 
well as the AGN strength dimension, we have plotted the average stellar population statistics for each of 
three groups (strong Type I, strong Type II, and weak Type I AGN) in Figure 20. Notice that the overall 
level of young and intermediate stellar populations is clearly different in the strong and weak populations, 
regardless of AGN Type. 

It is difficult to directly compare these observations with Kauffmann's finding that young stellar pop- 
ulations are correlated with AGN strength, as this result was obtained for both strong and weak Type II 
AGN. We do not sample weak Type II AGN in this paper. However, combining Type I and Type II popu- 
lations, there is a correlation between [OIII] line huninosity and the presence of dominating young and/or 
intermediate aged populations at any radii. This is consistent with Kauffmann's low-redshift finding that 
(1) AGN strength is correlated with the presence of young stellar populations provided that (2) the stellar 
populations of Type I AGN are indistinguishable from Type II AGN. Figure 21 plots [OIII] line luminosities 
for all AGN sources against the radially averaged fraction of young stellar populations (age < 100 Myr). 
This plot mixes Type I and Type II sources. 

The morphologies of the strong Type II AGN in this study follow trends outlined in Kauffmann's 
low-redshift sample. The morphologies of the galaxies in KaufFman's sample with considerable [OIII] (log 
Lioiii] > 7.8), a group of which this study's XID 155, 266, and 56 are members, fall into three broad cate- 
gories: (1) Single blue spheroids/amorphous galaxies, with no close companions and no significant structure. 
These look like E/SO galaxies with anomalous blue colors. (2) Single disk galaxies. In these, structure is seen 
outside the bright core, like arms, bars, and dust lanes. No close companions are visible in this category. 
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(3) Disturbed, interacting galaxies in which the core system possesses one or more close companions and 
obvious tidal debris. Of this study's strong Type II AGN, XID 56 falls into category (3), XID 266 falls into 
(2), and XID 155 falls into either (2) or (3), depending on whether the extended blue disk is a disk or tidal 
debris. As Kauffmann says, this shows that there is a variety of ways to feed the central black hole, but 
common to all of them are large patches of gas and extra star formation. 

We find that 3 of our 3 Type II AGN are LIRGs, as compared to 1 of 3 weak Type I AGN. 2 of 2 strong 
Type I AGN are LIRGS. In the Type I AGN sample (both strong and weak), the 24 micron to 8 micron flux 
ratios indicate that the source of the mid-IR emission is due to hot dust in close proximity to the central 
accretion disk itself. However, the strong Type I AGN are distinguished from the weak Type I AGN by the 
level of 24 micron flux, which elevates them to LIRG-status. The powerful Type II AGN, in sharp contrast, 
display 24/8 micron flux ratios indicative of a strong starburst contribution to the mid-IR emission, which 
is more than enough to classify them as LIRGs. It is unclear if the Type II AGN would remain LIRGs if the 
starbursting component were removed. These observations are consistent with the notion that starbursting 
in Type II AGN is somehow connected with their status as obscured AGN. All of the strong AGN (log 
L[oiii] > 7.5 in bolometric solar luminosities) are classified as LIRGS regardless of obscuration; however, 
certainly the unobscured AGN are LIRGs because of mid-IR emission associated with the AGN. 

There is additional information in the presence of point sources and their SEDs in the sample. The 
Type II population does not display point sources in any band, according to the criteria that the with 

a fitted point source be significantly (5(t) improved compared to that with a wider (100 mas) Gaussian. 2 
of 5 Type I AGN display clear point sources in all bands. These 2 AGN additionally display broad-lined 
emission. The SEDs of the point sources are consistent with those of unobscured QSOs. No other Type I 
AGN (XID 83) display point sources in any bands. 

Overall, the morphologies, stellar population distributions, and color distributions of the high-redshift 

sample of 8 AGN arc similar to the low-redshift samples, except for various small differences arguing for larger 
gas reservoirs in the strong Type II population. The evidence for this is drawn from two critical observations: 

(1) High-z strong Type II AGN are more likely to be LIRGS via starbursting events. This compares to all 
Type I AGN, which, if presenting LIRG activity, are excited to that level via hot dust surrounding the AGN. 

(2) Type II AGN hosts at high-z tend to be Sc/d/Irr galaxies, while Type I AGN tend to be earlier types 
(verified in Pierce et al 2007). With the findings that Type II AGN uniquely possess extended star formation, 
LIRG starbursting, and are morphologically distinct from Type I AGN, this raises the question: What is 
truly causing the obscuration that fundamentally distinguishes these two categories of AGN? If the obscuring 
medium is associated with these large scale starbursts, then this leads to the supposition that high-z Type 
II AGN are obscured via kpc-scale dusty features that are byproducts of starbursting in addition to small 
(r < 300 pc) dust tori, as in the case of Type I AGN. 

7.3. Significance of Stellar Populations Distinctions 

The discussion above must be treated with caution, as the number of galaxies in each bin (Strong Type 
I, Strong Type II, and Weak Type I) is either two or three. A student's t-test can be used to estimate the 
probability that the hosts in these groups have been selected from fundamentally different populations. We 
use the radially averaged fraction of young and intermediate aged stars (age < 1 Gyr), without weighting 
by luminosity, to compare the different populations. Comparing the Strong Type I population to the Weak 
Type I population gives a probability of 83% that the hosts are drawn from a different sample. Using just 
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the radially averaged fraction of young stars (age < 100 Myr) gives a probability of 91% that these two 
populations are distinct. However, when comparing the Strong Type I population and the Strong Type II 
population for stars less than 1 Gyr in age, the probability is 61%. This places higher weight on the conclusion 
that AGN strength {L^oiii]) is correlated with the presence of YSPs, but less weight on the observation that 
the stellar populations of Type I and Type II hosts at equivalent AGN strength are distinct. A larger, 
carefully-controlled sample of Type I and Type II AGN hosts at high redshift will be necessary to examine 
this question in more depth. 



8. SUMMARY 

We have completed a pilot study of a small sample of X-ray selected AGN in the GOODS field, combining 
K' band imagery from the Keck LGSAO system with deep ISAAC Ks imagery and HST ACS imagery in 
four optical bands. We have combined evidence from multiple independent AGN selection methods to verify 
that the eight galaxies in this sample are AGN. A morphological analysis of the galaxies compared to a 
sample of nonactive systems has revealed that the Type I AGN hosts are more likely to be ellipticals and the 
Type II AGN hosts are more likely to be spirals, disks, or irregular types. A stellar populations analysis of 
the five band photometry at a resolution better than 100 milliarcseconds in all bands has demonstrated that 
the stellar populations of the AGN hosts are consistent with their morphological identities, i.e., early types 
have older ages and later types have younger ages. We combine Type I and Type II populations and find 
that there is a correlation between AGN strength (^[o///]) ^-nd the presence of d ominating young and/o r 
intermediate aged populations at any radii, similar to what is seen at low-redshift ( Kauffmann et al. 2003h . 



The mid-infrared properties of the AGN hosts identify all of the strong Type II AGN hosts as LIRGs via 
star bursting events, the strong Type I AGN hosts as LIRGS via hot, localized dust surrounding the central 
engine, and the weak AGN as non-LIRGs. The strong Type II AGN hosts are also more likely to possess 
extended star formation than the Type I AGN. Although this study cannot come to rigorous conclusions 
about the nature of Type II AGN obscuration due to the small sample size, these last two pieces of evidence 
taken together with the previously seen correlation between morphology and AGN Type at z 1 suggest 
that Type II obscuration is due to kpc-scale dusty features that are associated with LIRG-level starbursting. 
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A. APPENDIX A: LGSAO PSF ESTIMATION 

In Laser Guide Star AO systems, the spatially and temporally variant PSF is more difficult to model 
than in an NGS system, due to additional error terms and the shifting atmospheric geometries (i.e., the cone 
effect). In a typical CATS field with LGS at the center and a tip/tilt star near the edge of the field, the PSF 
distribution is approximately given by the sum of three error components: (1) Anisoplanaticism and cone 
effect from the LGS itself, (2) Tilt anisoplanaticism (or anisokineticism), and (3) non-spatially dependent 
error sources such as shot noise from the tip/tilt star. LGS anisoplanaticism is assumed to be a radially 
symmetric distribution centered on the LGS. Anisokineticism is also radially symmetric, but is centered on 
the tip/tilt star. Error due to source (3) arc assumed to not depend on geometry. 

We use the Steinbring et al. (2005) formulation to estimate the PSF on a target-by-target basis using 
an exposure of a star cluster and a single stellar source in the target field. This semiempirical method 
assumes that the spatial and temporal components of PSF variation are largely separable. It utilizes an 
LGS image of a cluster taken once throughout the night to capture the spatial variance of the PSF (i.e., 
LGS anisoplanaticism and anisokineticism). Ignoring LGS dimming due to intervening clouds or technical 
problems with the laser, the remaining temporal variance throughout the night can be captured with an 
additional convolved spatially-variant Gaussian kernel to model image degradation due to anisokineticism 
from the tip/tilt. In the images presented, the tip/tilt star is imaged in the corner of the field to serve as 
the reference for the Gaussian. The tip/tilt star and laser guide star are not necessarily oriented in the same 
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geometry for cluster calibration images and science frames. Practically, this requires an estimation of the 
anisokineticism in the cluster image. 



In the lSteinbring et al.l (j2005[ ) method, PSFs are computed for each science target. A star is first chosen 
in the cluster image for which the orientation relative to the laser guide star is the same as the tip/tilt 
star in the science frame of interest. A second star is chosen in the cluster image close to the location of 
the particuliar science target relative to the laser guide star. The second star is deconvolved from the first 
star using the Max_Likelihood deconvolution function in IDL, which provides very good agreement upon 
reconvolution (< 5% at all pixels in a radial profile). The resulting kernel, which to first order translates 
PSFs at the tip/tilt star in the science frame to PSFs at arbitrary locations, includes the effects of LGS 
anisoplanatism and the difference in anisokineticism that these two cluster stars have seen. The tip/tilt star 
in the science frame is convolved with this kernel and further convolved with a Gaussian of full-width at 
half-max: 



Here Og is the angle between the science target and the tip/tilt star. Oi is the angle between the 
first cluster star (chosen to correspond to the science tip/tilt star) and the cluster tip/tilt star. 82 is the 
corresponding angle between the second cluster star (chosen to correspond to the science target) and the 
cluster tip/tilt star. 9a is the anisokinetic angle and R is a constant. The ratio {R/6a) is calibrated by a 
second star in the science frame. Note that the anisokinetic angle is assumed to b e the same between the 
cluster frame and the science frame. This is certainly not guaranteed to be true. Steinbring et al. ( 20051 ) 



recommends that a second calibration image of the cluster be taken with tip/tilt correction off (to obtain an 
independent measure of the anisokinetic angle), but these calibration images were not taken for the frames 
presented here. 

This uncertainty in the anisokinetic angle, which may vary by a factor of two, causes an uncertainty in 
the final Gaussian FWHM. We translate this into a quantitative PSF error via Monte Carlo simulation of 
multiple PSFs. Targets near the tip/tilt star will have less of this error because their error is dominated by 
LGS anisoplanatism instead. Additionally, this error is concentrated in the central 300 mas aperture because 
the Gaussian error functions due to anisokineticism are not widely extended. 

Steinbring et al. ( 20051 ) reports that the full method (with a both cluster calibration frames, with and 



without tip/tilt correction) gives 20% one-sigma error on the Strehl ratio on stars 8'.' from the LGS and 28" 
from the tip/tilt star. This is likely dominated by the uncertainty in the anisokineticism that we see. 



B. APPENDIX B: MORPHOLOGICAL PROPERTIES FOR COMPARISON SAMPLE 

Image tiles, GALFIT fits, model residuals, and morphological analyses are displayed for the normal, 
comparison sample in Figures 22-30. 



B.l. Normal 1 



Normal 1 is an elliptical galaxy, with smooth profiles in all bands. The only residuals remaining after 
subtraction of a best-fit sersic + Gaussian model are associated with the improperly-fit width of the central 
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core. The axis ratio of the sersic component is greater than 0.9 in all bands. There are elongated features 
at r = I'.'O in the optical residuals that resemble shells and rings. 

B.2. Normal 2 

Normal 2 is an inclined disk with a central bluish, elongated feature. As seen in the GALFIT residuals, 

the central structure possesses knot-like signatures of star formation. The morphology outside of the central 
region (r > 0'.'2) is smooth and disk-like in all bands, with no knot-like star formation signatures. This 
suggests that Normal 2 is an Sa or SO type galaxy. 

B.3. Normal 3 

Normal 3 is a nearly edge-on disk with a redder outer disk and indications of central star formation. 

The outer disk (r > 0'.'5 deprojected) is more prominent in bands redder than V. There arc two major central 
features, including a star forming ring of deprojected radius r = Ol'S and a central point source. The point 
source is visible in bands redder than V and the star forming ring is visible in bands bluer than z. 

B.4. Normal 4 

Normal 4 is an elliptical galaxy with smooth profiles in all bands. There is a strong central blue 
component, while the redder bands possess more extended profiles. The axis ratios of the primary sersic 
component in the best-fit sersic -|- Gaussian models are greater than 0.9 in all bands. There are few features 
in the residuals after subtracting galaxy models. 

B.5. Normal 5 

Normal 5 is a nearly face-on spiral galaxy (Sc) with clear evidence of star formation. Multiple blue 
star-forming knots are pinned to several spiral arms and a central star forming ring of radius 0'.'4. The 
central core is apparent in all bands, but all other knots and features disappear in the observed K', which 
only consists of the core embedded in an extended, spheroidal component. This spheroidal component is not 
apparent in the bluer bands. The two major spiral arms appear to emerge from the central core, but are not 
symmetric. 

B.6. Normal 6 

Normal 6 is a slightly inclined disk with a star forming ring and a central core. The central bulge is not 
visible in the B band, but is seen in all redder bands. The ring is composed of several pinned star-forming 
knots that are prominent in the bluer bands. Underlying the blue ring is a red, disk-like smooth structure 
seen primarily in the z and K' bands. The observed K' image is smooth and free of knots and rings, except 
for the central core. 
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B.7. Normal 7 

Normal 7 is a flattened elliptical (axis ratio ~ 0.7 in all bands) with a smooth morphology. 

B.8. Normal 8 

Normal 8 is a barred, grand design spiral galaxy with a central red bulge. The K' band (rest-frame 
J-band) image displays a spheroidal bulge on top of an elongated bar seen in the bluer bands as well, but 
does not show the spiral arms. The spiral arms are prominent in the bluer bands. The morphology is largely 
smooth, with an absence of star-forming knots and rings. 

B. 9. Normal 9 

Normal 9 is a compact, inclined disk with several star-forming knots in the central region (r < 0'.'5). 
The observed K' morphology is different from the optical bands, with two weak spiral arms embedded within 
an ellipsoidal structure. 

C. APPENDIX C: STELLAR POPULATIONS PROPERTIES FOR THE COMPARISON 

SAMPLE 

Stellar population fraction plots, radial profiles, and short discussions for the normal, comparison sample 
galaxies, as in Figures 12-19 for the AGN sample, are displayed in Figures 31-39. 

C. l. Normal 1 

This elliptical nonactive galaxy is dominated by older stellar populations at all radii. The shells and 
features apparent in the model-subtracted residuals, which were pointed out in §6, appear also to be composed 
of older stellar populations, meaning that they remnants of some long ago event. 

C.2. Normal 2 

Morphologically, this disk-like galaxy is an Sa/SO with a central blue core. The central core is composed 
of young stellar populations primarily. The region between r = C'S and r = 0'.'9 is a mixture of intermediate- 
aged and older stellar populations, and the region beyond this radius is dominated by an older population. 
This system has a positive age gradient. 

C.3. Normal 3 

This inclined disk is a mixture of populations of various ages, but intermediate ages dominate at small 
and large radii. Although this disk only possesses knot-like structure and a ring in the central 0'.'5, it is clear 
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that some star formation had occurred more recent than 1 Gyr to construct the outer disk. 

C.4. Normal 4 

This elhptical galaxy with a blue core and smooth features is composed of intermediate-aged populations 
in the central regions (r < 0'/4) and older stellar populations beyond this radius. 

C.5. Normal 5 

Normal 5 is a clear facc-on spiral with star-forming knots across the length of two spiral arms. Young 
stellar populations are dominant in the regions filled with the central star-forming ring, or r < 0'.'6. The spiral 
arms themselves, which emerge from the star forming ring at this radius, are described as an intermediate- 
aged component at this point, although the underlying older stellar population dominates at this and higher 
radii. 



C.6. Normal 6 

This irregular galaxy appears to have a red, bulge-like component in the center and a blue ring at 
r = (X'S — 0'.'4. The red core is largely an older stellar population, and the star forming ring clearly shows as 
a younger contribution. At higher radii, the system is a combination of equal fractions of populations, but 
the overlap of these points compared to the errors bars suggest that the photometric errors are too large to 
draw conclusions. 



C.7. Normal 7 

This smooth elliptical with very little blue contribution is dominated by older stellar populations at 
low and high radii (r < 0'.'3 and r > 0'.'7). Intermediate-aged populations dominate at radii between these 
two regions, but since there are no morphological features at these radii that might correspond to these 
populations, this might be a result of photometry error in the B-band. 

C.8. Normal 8 

Normal 8 is a clearly barred spiral with no knot-like structure and a core bluer than the surrounding 
system. This system appears to be the sum of three morphologically orthogonal components that are 
composed of different populations: A blue, star-forming component in the central core, an old bar, and 
intermediate aged spiral arms. The central core is an equal mixture of these different components. At higher 
radii (r > 0'.'2), the red bar dominates the photometry with an older stellar population. The spiral arms 
appear at the highest radius with an intermediate-aged contribution of 10% by mass. 
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C.9. Normal 9 

This inclined disk with the appearances of central star formation turns out to be dominated by a young 
stellar population for r < 0'.'6. There is a significant old component in the core, possibly due to the peak of 
some red spheroidal, bulge-like component. Beyond 0'.'6, the population is dominated by a weak, underlying 
old stellar population. 



- 29 - 



Fig. 1. — SEDs of point sources subtracted from the center of XID 32 (left) and XID 15 ( right 
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Fig. 2. — Log-log radial profiles for XID 56 (K') images at various intermediate steps during the calculation 
of aperture corrections. The solid line is the measured radial profile for the raw K' image. The dotted-dashed 
line is the radial profile of the best GALFIT fit to this galaxy, including a convolution to the known PSF. 
The diamond points represent the annulus photometry directly measured on the raw data. The dashed line 
is the radial profile of the GALFIT model without convolution with the PSF, denoting the best estimate of 
the intrinsic light distribution. Aperture corrections are computed by comparing the GALFIT model with 
and without convolution by the PSF. The triangle points plot the aperture-corrected photometry. 
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Fig. 3. — Tiled magery, GALFIT models, and colormaps for XID 56. The log-scale tiled images are 3'.' 
across. The five bands are split column-wise in the order B, V, i, z, and K'. The bottom five panels are 
the raw, inverted images. The middle five are the corresponding GALFIT models with bulge and disk sersic 
components. Residuals are displayed in the top row. The B-i colormap for XID 56 is shown in the bottom 
left and the V-K' colormap is shown in the bottom right. All magnitudes are in the AB system. K' imaging 
from adaptive optics is used within the gray square of 2.0" size; deep ISAAC imaging in Ks is used outside 
of the square. The correction from Ks to K' is small 0.02 mag) and ignored for the colormaps. 
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Fig. 7. — Tiled imagery, GALFIT models, and colormaps for XID 536, as in Figure 3. 
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Fig. 10. — Tiled imagery, GALFIT models, and colormaps for XID 32, as in Figure 3. 
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Fig. 11. — Measured SEDs and best-fit Bruzual & Chariot spectra for two cases: (left) the SED from the 
fourth ring of XID 155 and (right) the SED from the fifth ring of XID 266. The plotted spectra represent fits 
to the SED when the age is fixed at different values and the dust extinction and tau parameter are allowed 
to float. The solid line in both cases is the best-fit model to the SED when all parameters are allowed to 
float. Note that ages below 1 Gyr are largely ruled out in the left panel, but not so in the right panel, 
where models between 500 Myr and 4 Gyr are within the photometric error. The SED in the left panel is 
representative of non-degenerate cases and the SED in the right panel is representative of degenerate cases. 




Fig. 12. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 56. The 
radial profile plots the surface brightness for each band (AB magnitudes per square arcsecond) as a function 
of radius. The fractions of all stars formed by mass is plotted in a log scale versus the radius for three 
populations: Young (age < 100 Myr), Intermediate (100 Myr < age < 1 Gyr), and Old (age > 1 Gyr). The 
young population is denoted by blue asterisks, the intermediate age population is denoted by green triangles, 
and the old population is denoted by red diamonds. 
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Fig. 13. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 155, with 
plots and symbols as in Figure 12. 
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Fig. 14. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 266, with 
plots and symbols as in Figure 12. 
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Fig. 15. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 83, with 
plots and symbols as in Figure 12. 
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Fig. 16. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 536, with 
plots and symbols as in Figure 12. 
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Fig. 17. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 594, with 
plots and symbols as in Figure 12. 
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Fig. 18. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 15, with 
plots and symbols as in Figure 12. 
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Fig. 19. — Radial profiles and fractions of different stellar populations plotted vs. radius for XID 32, with 
plots and symbols as in Figure 12. 
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Fig. 20. — Plots of mean stellar populations statistics for three AGN populations: Strong Type II (upper 
left), Strong Type I (upper right), and Weak Type I AGN (lower left). The x-axis is the radial location in 
arcseconds and the y-axis is the fraction of stars in a given stellar population. Colors denote different age 
ranges of stellar populations, as explained in §7. 
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Fig. 21. — [OIII] line luminosities for all AGN sources plotted against the radially averaged fraction of young 
stellar populations (age < 100 Myr), not weighted by luminosity. 
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Fig. 22. — Tiled imagery and GALFIT models for Normal 1, as in Figure 3. 
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Fig. 23. — Tiled imagery and GALFIT models for Normal 2, as in Figure 3. 
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Fig. 24. — Tiled imagery and GALFIT models for Normal 3, as in Figure 3. 
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Fig. 25. — Tiled imagery and GALFIT models for Normal 4, as in Figure 3. 
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Fig. 26. — Tiled imagery and GALFIT models for Normal 5, as in Figure 3. 
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Fig. 27. — Tiled imagery and GALFIT models for Normal 6, as in Figure 3. 




Fig. 28. — Tiled imagery and GALFIT models for Normal 7, as in Figure 3. 
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Fig. 29. — Tiled imagery and GALFIT models for Normal 8, as in Figure 3. 




Fig. 30. — Tiled imagery and GALFIT models for Normal 9, as in Figure 3. 
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Fig. 31. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 1, with 
plots and symbols as in Figure 12. 
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Fig. 32. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 2, with 
plots and symbols as in Figure 12. 
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Fig. 33. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 3, with 
plots and symbols as in Figure 12. 
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Fig. 34. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 4, with 
plots and symbols as in Figure 12. 
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Fig. 35. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 5, with 
plots and symbols as in Figure 12. 
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Fig. 36. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 6, with 
plots and symbols as in Figure 12. 
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Fig. 37. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 7, with 
plots and symbols as in Figure 12. 
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Fig. 38. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 8, with 
plots and symbols as in Figure 12. 
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Table 1: Tabic of fields imaged with Keck LGSAO and NGSAO. 



Field 


Tip-Tilt Star Coord. 


Offsef^ 






Date 


Sources in Image 


GS 40s 


3h32ml3.8« 


-27°42'13'.'9 


S 


15.0 


20.6 


2005 Sep 28 


XID 266,56,83 


GS 33 


3h32i°55.5« 


-27°51'6'.'4 


W 


15.4 


22.2 


2005 Sep 28 


XID 15, Norm5-6 


GS 68w 


3h32ml0.3« 


-27°43'7'.'2 


NW 


16.0 


21.1 


2005 Sep 27 


XID 155,536,594,Norml-4 


GS 25 


3h32m36.1>5 


-27°40'5'.'0 


NE 


16.9 


21.7 


2006 Dec 11 


XID 32 


GS 71 


3h32m3 1.9*5 


-27°41'5'.'6 


W 


16.1 


21.0 


2006 Dec 11 


Norm7-9 



"R-band magnitude of guide star 

'distance between guide star and center of science field in arcseconds 
'^Direction from guide star to center of science field 



Fig. 39. — Radial profiles and fractions of different stellar populations plotted vs. radius for Normal 9, with 
plots and symbols as in Figure 12. 
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Table 2: Table of critical parameters of our AGN and comparison samples. 



XID #^ 


Coordinates (J2000)^ 


Redshift'' 


z-type'' 


B (AB)'^ 


R (Vega)^ 


X-ray c' 


83 


3h32ml4.98'5 -27°42'24'.'9 


1.76 


Photo 


25.13 


23.08 


AGN I 


155 


3h32m7.98« -27°42'39'.'4 


0.55 


Spect 


23.61 


22.37 


AGN II 


266 


3h32ml3.86s -27°42'49'.'0 


0.73 


Spect 


23.70 


21.70 


AGN II 


56 


3h32ml3.24s -27°42'4T/1 


0.61 


Spect 


21.68 


20.18 


AGN II 


536 


3h32ml0.76^ -27°42'34'.'6 


0.42 


Spect 


21.99 


19.39 


galaxy 


594 


3h32m9.71s -27°42'48'.'2 


0.73 


Spect 


24.71 


21.52 


AGN I 


15 


3h32m52.88« -27°51'20'.'0 


1.23 


Spect 


23.51 


22.62 


AGN I 


32 


3h32m37.47s -27°40'0'.'3 


0.66 


Spect 


23.52 


22.36 


AGN I 


Norml 


3ii32ml0.38^ -27°43'0'.'9 


0.422 


Photo 


23.12 


22.1 




Norm2 


3h32m9.23« -27°42'51'.'l 


0.369 


Photo 


24.64 


21.4 




Norm3 


3h32m9.40« -27°42'36'.'l 


0.729 


Photo 


25.43 


23.1 




Norm4 


3h32m8.41s -27°42'3T.'3 


0.563 


Photo 


25.02 


22.4 




Norm5 


3ii32m52.87^ -27°51'14'.'74 


0.892 


Photo 


24.03 


22.48 




Norm6 


31132^55.01^ -27°50'5T.'64 


0.607 


Photo 


24.41 


22.65 




Norm? 


3h32m31.65« -27°41'14'.'28 


0.764 


Photo 


25.72 


22.25 




NormS 


3h32m32.99s -27°41'17'.'06 


1.096 


Photo 


24.46 


22.57 




Norm9 


3h32m32.90^ -27°41'23'.'93 


1.076 


Photo 


24.76 


23.52 





"XID numbers, coordinates, and ground-based R-band magnitudes (Vega) for AGN are from'Giacco ni et al.l ||2002| V 
Spectroscopic redshifts are from iSzokolv et al. 1 1 (2004) and photometric redshifts are from Zheng et al] | |2004) . 

'^B-band magnitudes are measured directly from the ACS F435W (B) images, using a 3'.' diameter aperture. 
X-ray class is in lSzokolv et al. 1 1 |2004) and lZheng et al.l (l2004f) . 
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Table 3: S ummary of X- r ay, IR , and spectroscopic measurements for the AGN sample. The X-ray class is as 
defined in lSzokolv et alj (j2004l ). based on total 0.5 - 10 ke V X-ray luminosity and hardness ratio. The total 
0.5 - 10 keV X-ray luminosities presented here, taken from Tozzi et al. ( 2000(1 . a re in units of 10^^ ergs s~^ . 
The hardness ratio is defined as HR = (H-S) / (H-^R) as in lSzokolv et all |2004l ) and lCiacconi et al.l |2002h . 
where H and S are the net instrument count rates in the hard (2-10 keV) and the soft (0.5-2 keV) bands , 
respectively. The neutral hydrogen column densities (in units of 10^^ cm~^) are also from Tozzi et al. ( 20061 ) 
and are deri ved from fits to X -ray spectra. The X-r ay to optical f l ux rat io (denoted X/0) is as defined 



originally in iMaccacaro et al. I |l988,) and modified in ISzokolv et al.l (j2004l ) , with a correction for internal 



absorption in the Type II sources. The 24 micron flux (denoted /„ii(24)) utilizes 24 micron fluxes published 
in catalog form on the Spitzer GOODS webpage (http://data.spitzer.caltech.edu/popular/goods/, Dickinson 
et al., in prepar a tion). The 24 micron-to-8 micron flux ratio (denoted as 24/8), defined as log{i' f ^ {24) /i/fi, (8)) 



Brand et all ^M), 

additionally utilizes 8 micron fluxes extracted from the publically available images 



at the same website. Logarithmic [OIII] line luminosities are in units of the bolometric solar luminosity 
(4.0 X 10^3 ergs s"^). 



XID # 


-^X, 0.5-10 


HR 


Nh 


X/0 


/.(24) 


24/8 


log Loiii 




(xlO*^ ergs s^-^) 




(xlO^i cm-2) 


(dex) 


(A^Jy) 


(dex) 




83 


44.3 


-0.23 


< 2.0 


-0.1 


336 ± 7 


-0.18 ± 0.12 


7.9'" ± 0.4 


155 


3.07 


0.16 


35.9 


-0.6 


522 ± 7 


0.30 ± 0.08 


8.0^^ 


266 


35 


1.00 


887 


-0.7 


333 ± 6 


0.30 ± 0.10 


7.8^ 


56 


26.2 


0.11 


16.2 


-0.1 


774 ± 9 


0.37 ± 0.08 


8.5^^ 


536 


0.3 


-0.25 


< 3.0 


-2.3 


< 90^^ 


-0.40 ± 0.15 


6.3"" ± 0.5 


594 


2.44 


-1.00 


< 1.0 


-1.1 


< 77'i 


-0.22 ± 0.18 


7.0'" ± 0.5, < 6.8^= 


15 


63 


-0.45 


2.6 


0.2 


140 ± 4 




8.0'" ± 0.5 


32 


11.1 


-0.49 


1.3 


-0.1 


159 ± 4 


0.15 


7.3^= ± 0.2 



"Computed from OIII line intensities in lNetzer et al. 

^Lqiii estimated from X-ray luminosity, using fit for Type I AGN shown in figure 2 of lNetzer et al ] l|2006h . with appreciable 
scatter 

'^Lqiii estimated from spectra in lSzokolv et al. I ||2004|) 

''f 1,(24) measured directly from MIPS images at the Spitzer GOODS website above 
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Table 4: Summary of stellar populations analysis for AGN and comparison samples. The dominant stellar 
populations are listed for two locations: The inner region (r < 0.4") and the outer region (r > 0.7"). 
"Young" populations refer to those with an age less than 100 Myr, "intermediate" age populations have an 
age between 100 Myr and 1 Gyr, and "old" populations are older than 1 Gyr. 



ID 


X-ray class 


Morphology 


Inner Population, r < 0.4" 


Outer Populat 


56 


AGN II 


Irr 


young 


intermediate 


155 


AGN II 


Disk 


young 


old 


266 


AGN II 


Sc 


old 


intermediate 


83 


AGN I 


E 


intermediate 


old 


536 


weak AGN I 


E merger 


old 


old 


594 


weak AGN I 


E 


intermediate 


old 


15 


BLAGN P 


Disk 


intermediate 


intermediate 


32 


weak BLAGN P 


E 


old 


old 


Norml 




E 


intermediate 


old 


Norm2 




so/a 


young 


old 


Norm3 




Sa 


intermediate 


intermediate 


Norm4 




E 


old 


old 


Norm5 




Sc 


young 


old 


Norm6 




Irr 


old 


equal mix 


Norm? 




E 


old 


old 


NormS 




SBb 


equal mix 


old 


Norm9 




Sd 


young 


old 



"Broad line AGN from optical spectroscopy from lSzokolv et al.l l|2004l '). 



